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I  INTRODUCTION 

The  flow  field  around  the  Trailblazer  I  body  (  a  sphere  with  a  short  cylindrical 
afterbody)  has  been  calculated  for  three  flight  conditions  -  18,  900  feet  per  second 
at  150,  000  feet,  ZO,  ZOO  ft.  per  second  at  ZOO,  000  feet,  and  ZO,  400  feet  per  second 
at  Z50,  000  feet. 

In  the  study  of  flow  fields  about  blunt  bodies  at  these  speeds  and  altitudes,  not 
only  must  the  classical  assumption  of  a  calorically  perfect  gas  be  abandoned  but 
account  must  be  taken  of  various  chemical  reactions  which  become  significant  at 
these  conditions.  Each  reaction  introduces  a  characteristic  time  which,  if  com¬ 
parable  with  the  time  of  travel  through  the  region  of  interest,  will  cause  deviations 
from  thermodynamic  equilibrium  and  consequently  some  elaborate  description 
of  the  thermodynamic  behavior  of  the  fluid  is  required.  In  the  case  of  the 
relatively  high  altitudes  (above  150,  000  feet)  and  the  speeds  of  current  interest 
it  has  been  established  (Reference  3)  that  insofar  as  the  gross  thermodynamic 
properties  (i.e.,  pressure,  temperature,  etc.)  are  concerned,  it  is  sufficient 
to  consider  only  the  five  components  N,  »  O  and  NO  and  further  that  the 

equilibrium  concentrations  of  these  are  reached  relatively  close  to  the  bow 
shock  and  thereafter  remain  fixed.  This  assumption  governs  the  calculations 


reported  herein. 
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11  METHOD  OF  CALCULATION 

The  method  of  calculation  used  at  GASL  for  blunt  nosed  bodies  of  revolution 
has  been  outlined  by  Vaglio-Laurin  (Ref,  1)  and  Va glio- JLaurin  and  Ferri  (Ref,  2),  In 
short,  the  procedure  is  to  guess  a  shock  shape  in  the  nose  region  and  then  to  develop 
the  flow  field  in  the  slightly  supersonic  region  by  the  characteristics  method.  These 
data  serve  as  input  data  for  a  calculation  of  the  purely  subsonic  region.  The  result 
of  this  second  calculation  is  a  body  shape  which  may  be  compared  with  the  true  shape 
and,  if  necessary,  the  assumed  shock  shape  changed.  In  fact,  it  has  developed  that 
for  most  round  nosed  bodies  a  single  guess  will  suffice.  The  calculation  in  the  purely 
supersonic  region  then  proceeds  to  generate  the  flow  field  and  shock  shape  downstream 
of  the  nose.  These  calculations  have  been  programmed  for  the  IBM  704,  There  is  no 
way  in  these  programs  to  generate  a  shock  in  a  nnnuniform  flow  and  so  the  calculation 
must  stop  at  corners  or  surfaces  causing  shocks  within  the  rotational  region  downstream 
of  the  bow  shock.  In  the  case  of  the  flow  field  around  the  Trailblazer  I  body,  there  are 
two  such  shocks,  one  on  the  body  at  roughly  2  1/4  nose  radii  from  the  nose  and  the 
second  at  the  closing  of  the  wake.  If  it  is  desired  to  continue  the  calculation  beyond 
such  a  point,  the  shock  shape  must  be  estimated  and  conditions  downstream  of  it  used 
as  input  data  for  the  continuation  of  the  original  program,  (These  programs  are  being 
extended  to  automatically  generate  such  secondary  shocks  and  thus  avoid  the  laborious 
procedure  presently  being  used). 

The  general  procedure  outlined  above  is  independent  of  the  thermodynamic 
behavior  of  the  gas  about  which  some  a  priori  assumption  must  be  made.  This  assumption 
amounts  to  deciding  which  of  the  myriad  possible  reactions  are  significant.  It  was 
originally  decided  to  treat  nonequilibrium  effects  as  perturbations  on  a  basic  flow  which 
would  be  either  one  in  thermodynamic  equilibrium  or  one  of  constant  chemical  compo¬ 
sition  along  streamlines.  This  basic  flow  field  could  be  done  independently  and  then  the 
linearized  equations  obtained  by  neglecting  squares  of  the  perturbations  solved  to  get 
first  order  departures  from  the  basic  field.  This  proposed  method  was  abandoned  after 
it  was  shown  in  Reference  3  that  for  the  cases  of  interest,  the  mass  fractions  of  the 
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major  constituents  did  not  vary  (along  streamlines)  from  their  initial  values  behind 
the  bow  shock.  Furthermore  the  flow  in  the  subsonic  or  slightly  supersonic  region 
was  shown  to  be  in  thermodynamic  equilibrium  which  made  available  for  use  an 
existing  program.  The  degree  of  ionization  on  the  other  hand  could  be  calculated 
from  the  reaction  N+O^C^+e",  this  to  be  taken  in  equilibrium  at  the  temperature 
and  density  determined  by  the  frozen  flow  field.  As  mentioned  above,  the  justi¬ 
fication  for  these  simplifications  results  from  a  detailed  calculation  of  the  rate 
processes  in  Reference  3  and  no  further  discussion  will  be  given  here. 

With  these  assumptions,  the  basic  characteristics  equations  were  pro¬ 
grammed  for  the  IBM  704*.  The  detailed  procedure  used  was  then  to  assign  a 
shock  shape  and  generate  (with  the  assumption  of  thermodynamic  equilibrium)  the 
flow  in  the  transonic  region.  These  data  then  served  as  input  to  the  frozen  flow 
calculation  which  kept  the  composition  constant  along  streamlines.  Since  additional 
streamlines  entered  the  calculation  as  the  bow  shock  was  developed,  it  was  necessary 
to  calculate  the  chemical  composition  at  each  new  point  on  the  shock.  The  significant 
reactions  were  taken  to  be  N2^2N,  O2  :^20and  N+O  ^  NO  and  accordingly,  at  each 
point  on  the  shock  the  equilibrium  concentrations  of  these  five  were  computed. 


*  (Very  capably  by  Mr,  Edward  Lieberman  of  CASE) 
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III  RESULTS 

The  shock  and  streamline  shapes  are  shown  in  Figs.  1,  2  and  3  for  150, 

200  and  250,000  feet  respectively.  A  cylinder  (with  radius  1/2  that  of  the  nose) 
has  been  assumed  to  extend  downstream  from  the  wake  shock.  This  addition  is 
suggested  by  photographs  of  actual  wakes  where  the  v/ake  shocks  begin  at  roughly 
this  distance  from  the  axis  and  furthermore  and  possibly  more  to  the  point  the 
characteristics  equations  require  special  treatment  if  they  are  to  be  used  near  the 
axis. 

The  streamlines  sketched  in  these  figures  show  that  the  air  entering 
the  steep  portion  of  the  bow  shock  extends  over  a  large  portion  of  the  shock  layer 
downstream  of  the  body.  This  has  been  noted  earlier  (Ref.  4)  for  the  case  of 
thermodynamic  equilibrium.  Not  unexpectedly  then  the  density  profiles  in  the 
wake  show  a  relatively  constant  density  over  the  inner  portion  and  a  rppid  rise  in 
the  outer  portion. 

The  distributions  across  the  shock  layer  of  temperature,  density,  pressure 
and  electron  concentration  are  given  in  Figs.  4,  5  and  6  for  various  streamwise 
stations.  At  stations  on  the  body  these  exhibit  the  rather  unusual  characteristic 
that  the  temperature  is  lowest  on  the  body*  This  is  a  result  of  the  assumed  con- 
stant  composition  which  does  not  allow  the  energy  of  dissociation  to  be  released  as 
the  flow  expands  around  the  body.  The  increasing  temperature  and  decreasing  con¬ 
centration  of  atomic  species  in  the  direction  normal  to  the  axis  causes  the  electron 
concentration  profiles  to  have  maxima  at  these  stations. 

At  stations  in  the  near  wake,  the  profiles  exhibit  discontinuities  correspond¬ 
ing  to  the  trailing  and  wake  shocks.  (It  should  be  pointed  out  that  these  shocks  have 
been  estimated  both  in  position  and  strength  and  are  the  more  in  error  the  further 
from  the  body).  It  is  seen  that  these  shocks  cause  order  of  magnitude  increases  in 
electron  concentration  and  therefore  cannot  be  neglected. 
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In  the  far  wake  the  discontinuities  due  to  these  shocks  has 
disappeared  and  the  various  parameters  vary  monotonically  from  body  to 
shock. 

IV  CONCLUSIONS 

The  flow  field  around  a  typical  reentry  body  has  been  calculated 
with  the  assumption  of  frozen  chemical  composition  along  streamlines.  The 
effect  of  additional  shocks  in  the  flow  field  has  been  calculated  and  leads  to 
order  of  magnitude  changes  in  electron  density  in  the  wake.  It  is  thus 
apparent  that  some  effort  should  be  made  to  improve  the  accuracy  of  this 
part  of  the  calculation. 
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ADDENDUM  TO  GASL  TECHNICAL  REPORT  NO.  Zil 
TABLE  I  -MASS  FRACTIONS 


;am  lini! 

No. 

1 

1 

^0 

6/2. 

150, 

000  Ft. 

1 

. 00014 

.  235 

.  1480 

.  6110 

.  0062 

2 

. 00031 

.  233 

.  0963 

.  6627 

.0093 

3 

. 00035 

.  232 

.  0871 

.6718 

.  0097 

4 

.  0016 

.  228 

.0185 

.  7375 

.  0158 

5 

.  0088 

.  215 

.  0031 

.  7476 

.  0269 

200, 

000  Ft. 

1 

. 00005 

.  236 

.  2200 

.  5410 

.0028 

2 

. 00008 

.  236 

.  1941 

.  5669 

.  0033 

3 

. 00008 

.  235 

.  1572 

.  6027 

.  0042 

4 

. 00013 

.  235 

.  0975 

.  6620 

.  0050 

5 

.  0146 

.  212 

.  0008 

.  7507 

.  0218 

250, 

000  Ft. 

1 

. 00001 

.  237 

.  2460 

.  5160 

.  0012 

2 

. 00002 

.  237 

.  2190 

.  5427 

.  0015 

3 

. 00003 

.237 

.  1888 

.  5720 

.  0020 

4 

. 00004 

.  236 

.  1347 

.  6259 

.0023 

5 

. 00064 

.2  34 

.  0090 

.7499 

.  0059 

Attached:  Revised  Figures  I,  2  and  3 


